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MHD Simulations of a Moving Subclump with Heat Conduction 
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ABSTRACT 

High resolution observations of cluster of galaxies by Chandra have revealed the existence of an X-ray emitting 
comet-like galaxy C153 in the core of cluster of galaxies A2125. The galaxy C153 moving fast in the cluster core 
has a distinct X-ray tail on one side, obviously due to ram pressure stripping, since the galaxy C153 crossed the 
central region of A2125. The X-ray emitting plasma in the tail is substantially cooler than the ambient plasma. 
We present results of two-dimensional magnetohydrodynamic simulations of the time evolution of a subclump like 
C153 moving in magnetized intergalactic matter. Anisotropic heat conduction is included. We found that the 
magnetic fields are essential for the existence of the cool X-ray tail, because in non-magnetized plasma the cooler 
subclump tail is heated up by isotropic heat conduction from the hot ambient plasma and does not form such a 
comet-like tail. 
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I. Introduction 

High spatial resolution observations of central region 
of cluster of galaxies by Chandra have revealed the exis- 
tence of a peculiar X-ray emitting comet-like structure. 

Wang et al. (2004) presented results of a deep 
Chandra observation together with extensive multi- 
wavelength data of large-scale hierarchical structure 
related with A2125. An interesting feature is a dis- 
tinct X-ray tail on one side of the fast moving {v ^ 
1500 km s""'^) galaxy C153, probably created by ram 
pressure stripping. They suggested that C153 crossed 
the central region of A2125 containing cD-like ellipti- 
cal galaxies quite recently. Since X-ray emission above 
1.5 keV is absent in this tail, this tail is substantially 
cooler (fcr < 1.5 keV) than ambient intergalactic plas- 
mas [kT ^ 3.2 keV). They estimated that the length 
of this trail is ~ 22"(^ 88kpc) and its average width is 
~ 4"(~ 16kpc). Additionally, an extended [O II] line 
emission toward the same direction has been detected. 

Since heat conduction in cluster of galaxies can be 
very efficient (e.g., Takahara & Ikeuchi 1977), heat con- 
duction plays a key role in the formation of cold fronts 
(e.g., Ettori & Fabian 2000; Vikhlinin et al. 2001) and 
thermal balance in the cluster core (cooling flow prob- 
lem). Asai et al. (2004) showed the dramatic effect of 
magnetic fields on heat conduction in cluster of galax- 
ies. They carried out magnetohydrodynamic (MHD) 
simulations of a subcluster moving in a magnetized in- 
tercluster plasma. They showed that the contact sur- 
face (a cold front) between the cool subcluster plasma 
and hot intercluster plasma is maintained because the 
heat conduction across the cold front is suppressed by 
magnetic fields wrapping the forehead of the moving 



subcluster. In non-magnetized plasma, however, the 
cold front disappears by heat conduction. Similarly, 
the cool X-ray tail embedded in hotter ambient plasma 
should subject to heating by thermal conduction. 

In this paper, we investigate the interaction be- 
tween a moving subclump and magnetized intergalactic 
plasma and explore necessary conditions for the exis- 
tence of the cool X-ray trail. 

II. Simulation Model 

We simulated the time evolution of a subclump in a 
frame comoving with the subclump. We solve the two- 
dimensional (2D) resistive MHD equations in a Carte- 
sian coordinate system {x, y). We use the specific heat 
ratio 7 = 5/3. We assume an anomalous resistivity. 
The resistivity sets in locally only when the drift ve- 
locity (wd = p) exceeds the critical velocity {vc)^ 
where j is current density. When Vd > fc, the resistiv- 
ity is = r]o{V(i/vc — 1)^, otherwise ?] = 0. We adopt 
770 = 0.01 and Vc = 3.0. We set an upper limit of 
the resistivity, T^max = 1-0. We assume heat conduc- 
tion along the magnetic field line. The classical Spitzer 
conductivity (Spitzer 1962) is assumed. When mag- 
netic fields exist, the conductivity along the field line 
is K « K|| — kqT^/^{ko = 5 X 10~'^ergss~-^ cm~^ K~^) 
and the conductivity across the field line is k± = 0, 
where T is temperature. 

The size of the computational box is 600 kpc x 
400 kpc and the number of grid points is {N^, Ny) = 
(1200,800). 

The units of length, velocity, density, pressure, tem- 
perature, and time in our simulations are ro = 25 kpc, 
vq = 500kms"\ pq = 6.5 x 10"^''' gcm~^, po = 
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4 X 10 ^^ergcm ^, To = 1.5keV, and to = ro/i'o = 

5 X 10^ yr, respectively. 




Fig. 1. — Initial density distribution (color scale) in 
the central region (375 kpc x 250 kpc) of simulation box. 
Solid lines show magnetic field lines and arrows show 
velocity vectors. 

Figure 1 shows the initial density distribution in the 
central region (375 kpc x 250 kpc) . Solid lines and ar- 
rows show magnetic field lines and velocity vectors. 
The subclump is assumed to be a spherical isother- 
mal low-temperature (kTin — 1.5 keV) plasma confined 
by the gravitational potential of the subclump. It 
is embedded in the less-dense (pout = 0.25 po), uni- 
form hot (fcTout = 3.0 keV) plasma. Here the sub- 
scripts "in" and "out" denote the values inside and 
outside the subclump, respectively. We assume that 
the density distribution of the subclump is given by 
the /3-model profile, pin = Pc[l + (''/''c)^]"'^'^^^, where 
r = {x^ + y^y^^ and (3 = 2/3. The maximum den- 
sity is pc — 5po = 3.2 X 10"^^gcm~^, and the core 
radius is Tc — 8.3 kpc. The subclump is initially in hy- 
drostatic equilibrium under the gravitational potential 
fixed throughout the simulation. 

Table 1 shows model parameters. Important param- 
eter is the plasma beta (/3o) defined as the ratio of the 
ambient gas pressure to the magnetic pressure. When 

= Pgas/Pmag = 100, the magnetic field strength is 
B ~ 3pG. 

We assume that the ambient plasma initially has a 
uniform speed with Mach number M — Vx/cgout — 1, 
where Cgout is the ambient sound speed. The Mach 
number with respect to the sound velocity inside the 
subcluster is M' = Ux/csin = V^- 

Model HC is a non-magnetic model with isotropic 
heat conduction. Models MCa, MCb, and MCc are 
models with moderate magnetic fields (/3o = 100) and 
anisotropic heat conduction. In models HC, MCb, and 
MCc, the Mach number is taken to be M = 2, and 
in model MCa, it is taken to be M = 1. The inclina- 
tion of magnetic fields from motion of the subclump is 
parameterized hy 9 = arccos[v • B/(w B)]. 

We use a modified Lax-Wendroff method with ar- 
tificial viscosity for MHD part, and the heat conduc- 
tion term in energy equation is solved by the implicit 



red and black successive overrelaxation method (see 
Yokoyama & Shibata 2001 for details). 

For boundary conditions, the left boundary at a; = 
—5 is taken to be a fixed boundary, and other bound- 
aries are taken to be free boundaries where waves can 
be transmitted. 

III. Numerical Results 

(a) Effects of Magnetic Fields on Cold X-ray 
Tail 

Let us investigate effects of magnetic fields on the 
existence of the cold X-ray tail. Figure 2 shows the 
results for models MCa, MCb, MCc, and HC from left 
to right, respectively. The top and bottom panels in 
Figure 2 show the distributions of temperature and X- 
ray intensity at the central region (375 kpc x 250 kpc) 
after t = 2 x 10* yr. X-ray intensity is visualized from 
simulation results as logarithm of ~ p^. 

In models with magnetic fields (MCa, MCb, and 
MCc) , the cold plasma inside the subclump is expelled 
backward due to ram pressure and forms a cool tail 
dividing into two branches. The cool tail survives the 
heat conduction from hotter ambient plasmas, because 
magnetic fields wrapping the subclump suppress the 
heat conduction across them. On the other hand, in a 
model without magnetic fields (model HC), cool plas- 
mas inside the subclump is heated by thermal con- 
duction and evaporates quickly. Thus cool tail is not 
formed in this model. 

In all models with magnetic fields, magnetic fields 
accumulating ahead of the subclump form a magnetic 
shield (e.g., Miniati et al. 1999) and their strength is 
enhanced several times that of the initial state. In ad- 
dition to the suppression of heat conduction, magnetic 
fields also prevent the Kelvin-Helmholtz instability in 
this region. 

A typical length of the X-ray tail formed in our sim- 
ulation is 100 — 200 kpc, and its width is 25 — 50 kpc at 
t — 2 X lO^yr. The tails are longer in models MCb and 
MCc (Mach number M = 2), than model MCa (Mach 
number M = 1), because the subclump subjects to the 
stronger ram pressure. 



Table 1. Models and parameters. Heat conduction is 
isotropic in model HC. In other models, heat conducts 
only along magnetic fields. 9 is an angle between the 
motion of the subclump and the magnetic field. 
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Fig. 2. — Results for models MCa, MCb, MCc, and HC from left to right, respectively. The top and bottom 
panels show distribution of temperature (top) and X-ray intensity (bottom) at the central region (375 kpc x 250 kpc) 
after < = 2 x 10^ yr, respectively. Solid curves show the magnetic field lines and arrows show the velocity vectors. 



(b) Effects of Magnetic Fields on Energy Con- 
version 

When a subclump moves in magnetized plasmas, 
magnetic fields can extract the kinetic energy of the 
plasma moving with the subclump because magnetic 
field lines are stretched by the plasma moving with the 
sublcump. The kinetic energy is converted to magnetic 
energy and thermal energy. We computed the energy 
conversion rates in order to study whether a motion of 
a subclump heats the ambient plasma. Figure 3 shows 
the time evolution of magnetic (solid line), kinetic (dot- 
ted line), and thermal (dot- dashed line) energies with 
respect to the initial kinetic energy, respectively. The 
left and right panels show results for model MCb (with 
magnetic fields) and HC (without magnetic fields). 

The left panel shows that magnetic energy increases 
only slightly, because magnetic fields are deformed in 
small area close to the subclump. On the other hand, 
thermal energy increases while kinetic energy decreases. 
That is, the kinetic energy of the subclump is converted 
into thermal energy through a shock heating. The right 
panel shows the similar behavior in model HC. 

The inefficiency of the energy conversion from ki- 
netic energy to magnetic energy is partly due to the 
free boundary condition. When the magnetic fields 
are fixed at the boundaries, magnetic fields will be 
deformed until magnetic energy is comparable to the 
gravitational energy. 

IV. Discussion & Summary 



We carried out 2D MHD simulations of a subclump 
through a magnetized ambient plasma. In 
a)) we showed that the magnetic field is essential 




for the existence of a cold X-ray tail of a subclump like 
the galaxy C153 observed in A2125. Heat conduction 
across the magnetic fields is suppressed by the mag- 
netic field wrapping the subclump. This mechanism is 





Fig. 3. — Time evolutions of magnetic (solid line), 
kinetic (dot-dashed line), and thermal (dotted line) en- 
ergies integrated for the whole simulation region, re- 
spectively. Left: Results for model MCb (with mag- 
netic fields). Right: Results for model HC (without 
magnetic fields). 

the same as that which enables the maintenance of cold 
fronts in cluster of galaxies (Asai et al. 2004). 

In the context of interstellar matter, similar works 
have been done (e.g., Jones et al. 1996; Miniati et 
al. 1999), although heat conduction and gravity are not 
included in their simulations. 

They investigated, through 2D MHD simulations, 
the interaction of a uniform magnetic field oblique to 
a moving interstellar cloud. Miniati et al. (1999) dis- 
cussed the conversion rate of kinetic energy to magnetic 
energy for several models. In contrast to their results, 
only a small fraction of the kinetic energy is conv erted 
to the magnetic energy in our models (see ^IHrb)| | . The 
difference comes mainly from the fact that a subclump 
in our model has lower density and moves faster than 
that in their model. 

Makishima et al. (2001) proposed a model of heating 
of cluster plasma through the motion of member galax- 
ies in magnetic fields. In local simulations we presented 
in this papar, the energy conversion rate of kinetic en- 
ergy of the moving subclump to the magnetic energy is 
small because magnetic fields can freely move at bound- 
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aries. In cluster of galaxies, magnetic fields may be an- 
chored to the cD galaxies. When subclumps move in 
such magnetosphere of the cD galaxy, magnetic fields 
will be stretched and twisted. Under such situation, 
the kinetic energy of the dark matter clump will be 
extracted through the magnetic interaction. The de- 
formed magnetic field lines may form current sheets, in 
which magnetic reconnection converts magnetic energy 
into thermal energy and kinetic energy. We expect ef- 
ficient heating through this process. This mechanism 
can be the heat source which compensates for the ra- 
diative cooling in cluster plasmas. Obviously, we have 
to carry out global MHD simulations to study this pro- 
cess. We would like to report the results of such simu- 
lations in subsequent papers (Asai et al. 2005 in prepa- 
ration). 
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